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30MI  EFFECTS  OF  A  LATER©  SI3TEM 


m  DILATATIdAL  WAVES 


Abstraot 

In  this  report,  tha  effect  of  the  orustal  nodal  on  tha  raristion 
of  tho  surface  notion  with  tho  angle  of  inoidanoo  and  tha  frequency  la 
aacanLnad  for  aereral  orustal  nodal*  •  Tha  study  was  oarriad  out  hy 
programing  tha  problan  for  tho  IBM  1620  and  7072  computer  ays  tana  using 
tha  natrix  fomulation  originally  suggaatad  by  Thorn  on  and  parfsotsd  by 
Haakoll  and  Dcrnan.  From  thaaa  program ,  tho  ratios  of  tho  diaplaoe- 
nsnta  at  tho  fra*  aurfao*  to  the  total  aaplitude  at  dapth  war*  computed 
for  aararal  orustal  nodal  a  in  rang  as  of  fraquanoy  and  angla  of  inoidanoo 
of  intarast  in  seismology.  Thaaa  ratios  are,  in  of foot,  trananiaaian 
coefficients* 

Six  orustal  nodala  haring  such  features  as  thin  low-relocity  aur¬ 
fao  a  layers,  low -velocity  layers  at  depth,  and  relatiroly  thick  and 
rolatirely  thin  total  thicknesses  vara  considered.  For  aaoh  nodal,  tha 
transnission  coefficients  vara  computed  for  fraquanoias  ranging  fron 
•02  ops  to  10*0  ope  in  stops  of  .02  ops,  and  for  angles  of  inoidmoe 
ranging  from  21  degrees  to  5 3  degrees  in  stops  of  li  degrees.  Haskell's 
nodal  was  included  in  those  calculations  in  order  to  obtain  a  check  on 
the  calculations.  A  further  check  vas  obtained  by  using  the  transmis¬ 
sion  coefficients  to  synthesis*  the  surf ao*  motion  due  to  an  incident 
vara  of  the  form  ( 1/7 7)  (sin  20  77~t)/t  and  oonparing  these  rallies  with 
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those  predicted  1)7  ray  theory*  Tha  agrssnnt  vaa  wry  good* 

At  a  raault  of  these  calculation^  tha  inportnoe  of  the  frequency 
dependant  oharaoter  of  tha  oxuatal  affaot  haa  ban  furthar  emphasised* 
It  haa  ban  ah  own  that  a  thin  low-velocity  surface  layer  oauaaa  a  largo 
variation  In  tha  transmission  o oaf fio into,  while  a  low-velocity  layar 
at  dapth  haa  littla  affaot ,  aapaelally  at  low  frequencies  •  Further, 
it  haa  baan  shown  that  tha  total  orostal  thloknaas  la  ona  of  tha  aoat 
algnlfloant  f  aotora  In  determining  tha  variation  with  frequency.  Tha 
tranaadaaion  ooafflointa  of  Relatively  thick  orostal  models  hara  a 
much  aora  rapid  frequanoy  variation  than  thoaa  of  tha  ral  ativaly  thin 
orust*  At  frequencies  lass  than  .2  ops,  it  la  vary  difficult  to  dis¬ 
tinguish  between  oruatal  aodala  of  vary  nearly  tha  a  an  thloknaas  by 
nans  of  tha  frequency  variation  of  tha  transmission  eoaffioiants 
walaas  tha  intaxnal  atruotnre  is  vary  different. 
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son  mors  or  a  laxkd  soim  os  dilatjtkxul  vavis 

1*  Introduction 

With  tho  advent  of  largo  network  a  of  calibrated  seismographs,  It  has 
be  cone  Increasingly  Important  to  deter  nine  the  effect  of  the  layering  In 
the  earth’s  crust  beneath  the  individual  stations  on  the  eurfaee  notion 
generated  by  selamio  waves  striking  the  base  of  the  crust.  This  report 
studies  sons  aspects  of  this  problem  by  calculating  the  particle  notion 
which  would  be  observed  at  the  surface  of  a  system  of  plane,  horizontal, 
non-attenuating  layers  as  a  result  of  a  continuous  sinusoidal  plane  waive 
striking  the  base  of  the  layered  system  at  oblique  angles  of  Inoidanoe. 
(Fig.  1)  The  Haskell -Thomson  matrix  method  is  used  to  obtain  the  steady- 
state  reepoase,  and  then  the  transient  response  for  a  simple  pulse  shape 
is  obtained  by  Fourier  integration. 

Although  the  relatively  recent  matrix  formulation  is  used  in  this 
report,  the  observations  and  basic  theory  which  indicated  that  the 
orustal  layering  oould  have  a  significant  effeot  on  the  surface  motion 
have  long  been  known.  In  1929,  I manure  observed  small  oeoillations  af¬ 
ter  the  first  arrival.  He  suggested  these  oscillations  oould  be  due  to 
the  layering  in  the  earth's  crust.  Then  in  1932  Susuki  observed  a  varia¬ 
tion  in  the  apparent  angle  of  inoidenoe  with  f  requenoy.  Since  these 
phenomena  oould  not  be  explained  on  the  basis  of  a  homogeneous  half  space, 
various  Japanese  investigators  began  theoretical  studies  to  determine  the 
effect  of  layering  in  the  orust  on  body  waves. 
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Under  the  conditio—  of  continuity  of  stroea  and  displacement  at  the 
layer  boundaries,  ouch  investigators  aa  Saaava  (1930),  Hasegawa  (1930), 
Sesava  and  Xanai  (1932a,  b,  e,  193bj  1935,  1936)  and  Miahi— re  and  Kanai 
(1933)  ahowed  that  layering  in  tha  crust  aould  oaosa  tha  affaeta  — n- 
tie— d  ab ora*  Although  in  principle  thia  faet  vaa  wall  aatabllahad  by 
1936,  in  practice  little  vaa  known  quantitatively  about  the  magnitude  af 
the  affect  af  tha  layering  m  tha  aurfaea  motion.  Tkcaa  oaloulatie— 
which  ware  carried  eut  uaually  involved  aueh  si— lifying  aasumptic—  aa 
j  ncr—1  inoidanoa,  equality  of  dmaity,  red— ties  of  tha  problem  to  two 

layera  (i.e.  —a  layer  over  a  half-spade),  and  a  value  of  Faience '• 

;  ratio  equal  to  0,25*  However,  from  theae  atudiaa,  oapaoially  thoae  of 

\ 

Soaawa  and  Kanai  (193b,  1935),  it  vaa  pooaible  to  atata  that  tha  effect 
of  the  layering  varied  with  both  fraquanoy  and  angle  of  incidence. 

In  1939  Hiehinura  and  Takuya—  (1939a,  1939b)  oo— idered  tha  off— t 
of  a  a  ingle  two-layer  aye  tan  for  nine  anglea  of  inoidanoa  with  thirty 
values  of  ateady-etate  fraquanoy  at  each  angle.  This  work  repree— tad 
tha  a— t  extensive  sat  of  calculations  in  exist— aa  until  1962, uhen  M. 

A.  Haskell  published  his  work.  In  this  interval  tha  observations  rf 
Outanberg  (1957)  and  Muttli  and  Whitmore  (1961)  gave  even  —re  eephaala 
to  tha  importance  of  the  crustal  affect  — d  its  variation  with  frequency, 
as  did  tha  —dal  atudiaa  of  Proas,  Oliver  and  Being  (195b),  Ivanova 
(1959)  and  Lavergna  (1961).  In  addition,  investigators  such  —  Lindsay 
(1939),  Kanai  (1953,  1957),  Kanai  and  Teshisawa  (1959)  oarxied  —  ether 
theoretical  studies.  However,  tha  theory  and  — thod  of  approach  re¬ 
mained  basically  tha  s— a,  and  the  primarily  qualitative  so— lusis— 
may  be  sun— vised  aa  folio— s 


I)  The  small  oe  dilations  after  tha  arrival  of  a  body  phase  aan 
ba  oaused  by  tha  orustal  structure  rathar  than  tha  source 
motion, 

t)  Tha  affaot  of  tha  eruat  varies  with 

a)  oruatal  atruotura 

b)  angle  of  inoidenoe 
e)  fraquanoy, 

3)  Klliptioal  part lo la  motions  from  dilatational  aeuroaa  aan  ba 
explained  on  tha  baa la  of  layering  of  tha  eruat* 

II)  Low-velocity  sedimentary  aurfaoa  layers  oan  eaaaa  large  ampli¬ 
tudes  to  ba  obs erred  at  tha  free  aurfaoa. 

5)  Tha  affect  of  tha  oruat  la  nearly  ooaatant  for  wares  with  wary 
long  ware  length. 

These  conclusions  point  out  tha  need  to  oonaidar  tha  affaot  of  tha  oruat 
on  any  solsXLo  study  uhioh  is  based  on  amplitudes  or  on  tha  aasooiatad 
apparent  angla  at  inoidencs,  especially  since  tha  obaarratians  are  often 
serried  out  on  instruments  with  different  fraquanoy  response  ourres* 

Soma  seismologists  hare  attested  to  use  tha  fraquanoy  dependence 
of  tha  surface  motion  on  a  layered  orust  as  a  naans  of  determining 
orustal  structure.  In  particular,  sereral  investigators  hare  atteapted 
to  use  tha  speotrun  of  tha  ground  notion  to  determine  tha  orustal  layer¬ 
ing  (Clewell  and  Simon,  1950,  Barton,  1956a,  bj  Barton  sad  Spinatyeva, 
1950)  and  Khudtinsky,  1961)*  Howarar,  baosusa  of  tha  laok  of  ooqparablm 
quant it at ire  theoratioal  studies,  their  work  is  mainly  of  interest  as 
ons  of  tha  first  attempts  to  apply  tha  fraqusaoy  rariatien  of  tha 
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surface  motion  as  a  tool  to  determine  crustal  structure. 

In  addition  to  the  lack  of  quantitative  theoretical  studies, the  above 
work  suffered  from  its  inability  to  determine  whether  the  frequency  vari¬ 
ation  of  the  spectrum  of  the  ground  notion  arises  fron  the  crustal  inter¬ 
ference  or  from  the  source  motion. 

As  will  be  seen  in  section  2,  the  effect  of  the  source  can  be  elim¬ 
inated,  or  at  least  minimized,  by  considering  the  ratio  of  the  horizontal 
surface  motion  to  the  vertical  surface  motion  as  a  function  of  frequency. 
Some  attempts  to  folios  this  line  of  invesgation  were  carried  out  both 
theoretically  and  experimentally  by  Malinovskaya  (1959),  Ivanova  (I960); 
Halperin  and  Folova  (1960)|  and  Halperin  (1962).  Although  this  method 
appeared  to  be  superior  to  the  use  of  the  speotrua  of  a  single  component, 
these  investigations  again  suffered  from  a  lack  of  quantitative  theoreti¬ 
cal  values  which  could  be  used  for  comparison. 

two  events  combined  to  provide  a  means  of  obtaining  theoretical  quan¬ 
titative  values  (xi  a  sufficient  scale  to  be  able  to  separate  the  effects 
of  the  model,  the  angle  of  incidence,  and  the  frequency.  These  weres 
a)  the  development  of  the  matrix  formulation  of  the  boundary  condition, 
and  b)  the  general  availability  of  high  speed  computers. 

The  present  state  of  the  matrix  formulation  of  the  boundary  conditions 
is  primarily  due  to  three  investigators:  Thomson  (1950);  Haskell  (1953); 
and  Dorman  (1962).  Thomson  made  the  original  formulation  in  which  he 
treated  the  problem  of  propagation  of  plane  waves  through  a  system  of 
plane  parallel  layers.  Haskell  (1953)  corrected  an  error  in  Thomson's 
work  and  emphasised  the  application  of  the  method  to  the  dispersion  of  Love 


and  Rayleigh  waves.  In  1962  Dorman  showed  that  the  previous  formulations 
had,  in  offset,  over-specified  the  boundary  conditions  at  a  boundary  be¬ 
tween  a  liquid  and  a  solid.  Although  the  applications  to  surface  waves 
received  the  major  attention  of  seismologists  until  I960  and  1962  when 
Haskell  considered  the  applications  to  body  waves,  the  original  formula¬ 
tion  was  quite  general.  As  early  as  1953,  Matumoto  considered  the  trans¬ 
mission  of  body  waves  using  matrix  techniques.  At  that  time,  even  with 
this  powerful  tool  he  was  not  able  to  oonsider  as  many  cases  as  had 
Nlshimura  and  Takayana  (1937,  1939a).  However,  the  matrix  techniques  are 
well  suited  to  modem  digital  conputers  and  with  their  coming  it  has 
become  possible  to  treat  systems  of  many  layers  for  many  angles  of  in¬ 
cidence  and  many  steady-state  frequencies. 

In  1962  Has  kail  obtained  contour  plots  of  the  amplitude  and  phase 
to  the  vertical  and  horisontal  surface  motions  for  a  two-layer  system 
for  angles  of  incidence  ranging  from  sero  to  90  degress  and  for  fre¬ 
quencies  ranging  from  .01  to  *333  ops.  In  1963  Phinney  and  Smith  ob¬ 
tained  theoretical  ratios  of  the  spectra  in  the  low  frequency  range  for 
a  number  of  crustal  models.  In  196U  Phinney  showed  a  series  of  curves 
illustrating  the  results  of  a  large  number  of  calculations  cf  the  theo¬ 
retical  ratios  of  the  vertical  surface  motion  to  the  horisontal  surface 
motion  as  a  function  of  frequency. 

These  calculations  are  undoubtedly  the  most  extensive  set  of  cal¬ 
culations  which  have  been  done  in  the  frequency  range  from  0.001  to 
0.2  ops.  From  these  calculations  he  oonoluiods 
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1)  The  effects  of  intermediate  and  deep  crustal  structure 
is  isolated  in  the  low  frequency  range* 

2)  The  positions  of  peaks  in  the  ratio  do  not  ohange  appre¬ 
ciably  with  changes  in  the  angle  of  incidence* 

3)  Fine  structure  in  the  transfer  ratio  is  due  to  the  structure 
of  the  upper  mantle* 

b)  Thin,  low  Telocity  surface  layers  hare  no  singular  effects 
on  the  crustal  peaks  in  the  low  frequency  range* 

He  compared  these  commutations  with  the  obserrations  of  the  ratios  of 
the  spectra  of  the  ground  action  at  Albuquerque,  Mew  Mexico,  and  Bermuda- 
Coluabia.  From  his  analysis  he  was  able  to  obtain  crustal-upper  auntie 
models  which  were  geologically  satisfying*  Thus,  in  addition  to  obtain¬ 
ing  considerable  information  about  the  theoretical  variation  of  the 
spectral  ratios,  he  was  able  to  obtain  and  use  observational  studies 
for  the  determination  of  crustal  structure* 
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Although  the  matrix  formulation  of  the  boundary  conditions  in  a 
layered  system  is  of  considerable  importance  at  the  present  time,  it 
vill  not  be  developed  in  detail  here  beoause  of  the  many  excellent  tre*t- 
nents  which  already  exist  (Haskell,  1953}  Dorman,  1963).  However,  sinoe 
the  application  of  this  formulation  to  body  waves  has  reoently  received 
more  attention,  we  will  briefly  consider  the  theoretical  background  at 
this  application. 

As  a  result  of  his  analysis  in  1953,  N.  A.  Haskell  obtained  the 
following  equation  governing  the  steady-state  transmission  of  plane  si¬ 
nusoidal  waves  through  layered  systems t 


where  A  and  (0  are  plane  wave  solutions  to  the  dilatational  and  rota¬ 
tional  wave  equations.  The  singly  primed  terms  are  the  solitudes  of 
the  downward  travelling  components  while  the  doubly  primed  terms  corres¬ 
pond  to  the  upward  travelling  waves. 

is  a  It  x  1)  matrix  associated  with  the  boundary  condition 
between  the  half -space  and  the  n-lto  layer.  This  matrix 
is  independent  of  frequency. 
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Ajj-l  is  a  li  x  li  matrix  associated  with  trans  mission  across  ths 
n-1^1  layer  and  the  boundary  between  the  n-1^  layer  and 
the  n-2nd  layer.  Matrices  of  this  type  depend  on  the  layer 
constants,  the  angle  of  incidence  and  the  frequency  of  the 
input  plane  wave* 

u0  Is  the  horizontal  conponent  of  the  partiole  Telocity  at 
the  zeroth  interface  (Fig.  1) 

$0  is  the  vertical  conponent  of  the  partiole  velocity  at  the 
seroth  interface. 

c  is  the  apparent  surface  veloaLty. 

is  the  normal  stress  at  the  zeroth  interface, 
is  the  tangential  stress  at  the  zeroth  interface. 

3n  is  the  dilatational  wave  velocity  in  the  nth  layer. 

This  equation  relates  the  motions  at  the  top  and  bottom  of  a  layered  sys¬ 
tem,  incorporating  the  conditions  of  continuity  in  stress  and  velocity  at 
each  interface.  In  the  particular  case  being  considered,  there  is  no 
up-coming  shear  wave  in  the  nth  layer,  and  the  top  of  the  layered  system 
is  boiaaded  by  a  free  surfaoe.  Following  Haskell  (1962),  we  then  eet 

U)*  ■  0 
n 

CT0  -  o 

1  o  ■  0 

corresponding  to  the  above  conditions. 

Under  these  conditions,  we  can  solve  the  four  algebraic  equations 
which  are  implied  in  the  single  matrix  equation  for  An,  UJn,  u^/o, 
in  terns  of  An.  In  particular,  we  can  obtain 
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2(J32  -  J^)  fan 


2  *Jlil  ’  J31* 


where  the  J»s  are  the  elements  of  the  U  x  U  matrix  defined  by 
Kn“1  An_1  .  •  •  and  Den  -  (Jn"J21^J32"JU2^  “  ^12*^22 
Rather  than  use  the  ratios  of  Uq  and  $0  to  o  «e  developed  in  equa¬ 
tions  (2),  we  shall  find  it  more  convenient  to  obtain  the  ratio  of  the 
components  of  the  surface  particle  velocity  to  the  total  particle  veloc¬ 
ity  in  the  nth  layer  due  to  the  incident  dilatational  wave* 


A2- 


°LA 


Therefore  from  equations  (2)  and  (3)  we  have 


2  fa* 

«n 


<w 


These  same  equations  represent  the  ratios  of  the  surface  displacements 
to  the  displacements  at  depth  since  we  are  considering  steady-state 
si^>le  harmonic  notion.  Thus  we  have 
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Ug/ATn  ■  (J|t2-J32)/D«a 


(5) 


where  ATn  is  the  amplitude  of  the  total  particle  displacement  of  the  in- 

i.  L 

cident  dilatational  wave  in  the  n  layer*  Equations  (5)  ere,  in 
effect,  the  equations  for  the  transmission  coefficients  which  relate  the 
motion  at  the  free  surface  of  the  layered  system  to  the  motion  at  the 
base  of  the  system  directly  beneath  the  point  on  the  surface*  These  co¬ 
efficients  are  complex  functions  involving  the  frequency,  the  horizontal 
phase  velocity  (ihioh  Is  related  to  the  angle  of  incidence),  and  the 
layer  parameters.  We  shall  use  the  symbols  TU  and  TW  as  defined  below 
for  these  coefficients* 


(6) 


In  a  later  section  we  shall  examine  the  values  of  TU  and  TW  which  were 


obtained  by  evaluating  equations  (5)  for  a  mufcer  of  crustal  models* 
These  models  have  such  features  of  seismological  interest  as  relatively 


thin  low  velocity  layers  both  at  the  surface  and  at  some  depth,  and 
relatively  thick  and  relatively  thin  total  crustal  thicknesses*  Before 
doing  so,  however,  we  will  present  some  applications  of  TU  and  TW 


which  are  of  interest* 
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Since  TU  and  TW  are  dependent  on  cruatal  atrueture,  it  may  be 
possible  to  uee  maasureaents  of  thoee  quantitiee  to  determine  cruatal 
■tructure.  If  the  input  ware  is  of  the  fora  (AfB)  thma  the 

eurfaoe  notion  ean  be  repreaented  aa 


»o  "  U*n)  (T«)  ei(*>t-kx) 
«o  “  <*n)  (TW)  ei(p^"kx) 


Since  every  point  on  the  free  surf ace  is  executing  the  same  type  of 
notion,  although  with  different  phases,  we  oan  consider  the  point  x  -  o 
without  any  loes  of  generality*  Kith  this  understanding,  and  sinoe  the 
system  is  linear  we  oan  express  the  surfaoe  motion  due  to  an  incident 
wave  of  the  form 
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where  Jffn  is  the  frequency  spectrum  of  the  input  wave,  and  nay  be  a  com¬ 
plex  function  of  frequency. 

In  actual  practice  u^t)  and  w0(t)  are  recorded  by  seismographs  and 
from  this  notion  we  would  like  to  obtain  information  about  the  crust, 
or  possibly  about  the  aouroe.  Upon  examining  equations  (8)  we  see  that 
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the  ef foots  of  the  cruet  and  source  art  interconnected  in  the  quantities 
we  usually  observe.  In  order  to  separate  these  two  effects ,  let  us  con¬ 
sider  the  Fourier  transforms  of  the  ground  notion*  If  we  define  the 
transform  of  £  as  F(£),  we  have 

nu0(t))  .  (jib)  (to)  (» 

F(» 0(t))  •  («n)  (IW) 

The  left-hand  side  of  equations  (9)  can  be  determined  froai  the  observa¬ 
tions.  Thus  we  can  determine  the  product  of  the  souroe  spectrum  and  the 
transmission  effect  of  the  layers.  However,  if  we  take  the  ratio 
F(u0(t))/F(w0(t))  we  obtain  a  quantity  shioh  is  independent  of  the  fre¬ 
quency  content  of  the  source. 

rK  (O)  "  JU_  (10) 

F(w0(t))  W 

This  ratio  depends  on  the  frequency,  the  horizontal  phase  velooity 
(angle  of  incidence),  and  the  layer  parameters.  We  might  look  upon 
TU/TW  as  the  tangent  of  an  apparent  angle  of  incidence  Uhloh  is  a  func¬ 
tion  of  frequency. 

Since  this  ratio  depends  on  the  etructure,  it  should  be  possible 
to  determine  the  etructure,  or  the  clase  of  structures,  tfcich  give 
theoretical  ratios  to  matoh  the  ratios  obtained  from  the  observations 
at  a  given  site. 

Phinney  (1961*)  studied  the  crust  near  Albuquerque,  M.  M.,  by 
matching  ourves  by  trial  mad  error,  and  found  that  structures  which 
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vtn  in  agreement  with  the  known  geology  gar*  good  results,  while  those 
whioh  were  radio ally  different  gave  poor  agreement.  Thus  even  at  this 
tine,  this  method  of  determining  crustal  structure  appears  to  hold 
some  promise.  In  the  future  it  is  quite  possible  that  co^iuter  inver¬ 
sion  programs  similar  to  those  used  in  aurfaoe  wave  work  may  be  used. 

This  method  has  the  advantage  of  allowing  the  determination  of 
crustal  structure  beneath  a  station  from  the  observations  at  that  single 
station  alone.  Further,  it  does  not  depend  on  the  pioking  of  a  single 
event,  but  uses  the  series  of  events  whioh  are  reoorded  at  the  surface 
due  to  the  direct  arrival  of  the  plane  wave  plus  all  of  the  arrivals 
coming  from  crustal  reflections  and  conversions.  This  use  of  time  sec¬ 
tions  rather  than  specific  events  may  be  of  considerable  help  in  these 
esses  where  it  is  difficult  to  identify  specific  events. 

Once  the  crustal  structure  beneath  a  station  is  known,  then  it 
should  be  possible  to  determine  the  frequency  content  of  the  wave  which 
struok  the  base  of  the  crust.  Referring  to  equations  (9),  vs  see  that 


F(Uo(t))  .  F(w0(t)) 

TU  TW 


(11) 


This  would  allow  us  to  obtain  a  better  description  of  the  source.  A 
technique  such  as  this  may  prove  to  be  especially  valuable  in  shear 
wave  inveatigations  • 

It  should  be  emphasised  that  in  these  considerations  there  are 
some  important  restrictions  whioh  may  severely  limit  the  application 
of  these  techniques.  For  example,  the  attenuation  of  high  frequencies 
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which  will  pro ride  the  greatest  resolution,  the  crust  et  the  souroe, 
the  curvature  of  the  reel  earth,  and  the  necessity  of  having  perfectly 
Hatched  instruments,  all  have  effects  which  are  not  considered  here* 

At  the  present  tine,  we  oan  only  point  to  the  work  of  Phinney  (196k) 
which  indicates  that  these  aethods  are  applicable  at  least  to  the  study 
of  the  ooarse  features  of  orustal  layering* 

In  the  remainder  of  this  report,  we  will  examine  the  theoretical 
variation  of  the  transaission  coefficients  TU  and  TV  with  frequency, 
angle  of  incidence,  and  orustal  model  in  order  to  obtain  sobs  under¬ 
standing  of  the  offsets  of  these  pareastere* 
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3. .  Steady-State  Transalaaion  Coefficients 

In  order  to  study  ths  variation  of  TU  and  TV  with  frequency,  angle 
of  incidence,  and  crustal  structure,  wa  wrote  program  tar  tha  IBM  1620 
and  7072  computer  systems.  The  basic  program  in  1620  form  is  shewn  in 
Appendix  I  (the  7072  fora  of  this  program  differs  only  in  the  input - 
output  statements  and  the  addition  of  a  Fourier  integral  routine)*  The 
input  to  this  program  specifies  the  ranges  of  angle  of  incidence  and 
frequency,  and  the  thicknesses,  velocities  and  densities  of  tha  layers 
of  the  crustal  nodal  to  be  conn ida red.  In  this  study,  the  angle  of 
incidence  was  varied  froa  21  degrees  to  53  degrees  in  steps  of  k  de¬ 
grees,  and  the  frequency  was  rariad  from  *02  ops  to  10*0  ops  in  steps 
of  *02  ope*  These  ranges  were  chosen  bessusa  of  their  appl inability 
in  various  seisnologioal  studies* 

The  models  Wiioh  were  chosen  are  listed  in  Table  I  together  with 
a  brief  description  of  their  predominant  features*  Us  determined  the 
values  at  the  amplitude  and  phase  of  TU  and  TV  for  each  of  these  models 
at  the  values  of  the  angle  of  incidence  and  frequency  given  above*  Due 
to  the  extent  of  the  data  (greater  than  27,000  individual  points), 
only  a  representative  number  of  curves  will  be  shown  to  illustrate  the 
conclusions*  The  complete  set  of  data  is  available  in  tbs  form  of 
punched  cards  at  the  St.  Louis  University  Institute  of  Technology, 
Department  of  Geophysics.  It  should  be  noted,  however,  that  unless 
the  entire  range  of  data  is  desired,  it  is  often  more  convenient  to 
ooMpute  the  data  in  the  desired  range  by  using  the  program  in 
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Model  Description  of  Thickness 

P- 

s- 

Density 

Mo. 

Main  Feature  (Km.) 

Velocity 

(Km/Sec) 

Velocit; 

(Km/Sec 

I  (Ora/Cm^ 

1? 

A  crustal  model  1.0 

4.40 

2.50 

2.70 

similar  to  that  19.0 

6.20 

3.50 

2.80 

of  the  central  18.0 

6.40 

3.70 

2.90 

U.S.  with  a 
thin,  low- 
velocity  sur¬ 
face  layer 

8.20 

4.60 

3.30 

• 

2. 

The  same  model  20.0 

6.20 

3.50 

2.80 

as  1.  but  the  18.0 

6.40 

3.70 

2.90 

surface  layer 
is  removed  and 
the  second  layer 
is  extended  to 
keep  the  total 
crustal  thickness 
constant 

8.20 

4.60 

3.30 

3. 

A  crustal  model  15.0 

6.20 

3.50 

2.80 

with  a  relatively  5»° 

6.00 

3.30 

2.80 

thin  low- velocity  18.0 

6.40 

3.70 

2.90 

layer  at  depth 

8.20 

4.60 

3.30 

4. 

The  model  of  37.0 

6.285 

3.635 

2.869 

Haskell  (1962) 
included  as  a 
check  on  calcu¬ 
lations 

7.960 

4.600 

3.370 

5. 

A  crustal  model  30.0 

6.20 

3.50 

2.80 

representing  a  20.0 

6.40 

3.70 

2.90 

relatively  thick 
crust,  having  the 
same  velocities 
as  2  sind  6 

8.20 

4.60 

3.30 

5. 

A  crustal  model  15*0 

6.20 

3.50 

2.80 

representing  a  5.0 

6.40 

3.70 

2.90 

relatively  thin 
crust,  having 
the  same  veloci¬ 
ties  as  2  and  5 

8.20 

4 

4.60 

3.30 

Models  1  and  2  were  suggested  by  Mr.  T.V,  Mofivilly  as 
corresponding  to  crustal  structure  In  Missouri  and/or 
the  central  United  States. 
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Appendix  I, 

In  discussing  the  figures  on  the  following  pages,  we  shell  neke 
use  of  the  tern  "character"  and  "shift,"  "Character"  Is  Intended  to 
refer  to  the  relations  between  the  asplitudea  for  Yarlous  frequencies 
for  a  given  angle  of  inoidenoe  and  a  given  model.  That  is,  "charac¬ 
ter"  refers  to  the  general  shape  of  adjacent  peaks  and  troughs*  On  the 
other  hand,  "shift"  is  intended  to  refer  to  the  fact  that  the  fre¬ 
quency  variation  often  appears  to  be  c oppressed  or  extended  and  thus 
segments  of  the  curves  having  similar  character  are  displaced  in  the 
frequency  domain* 

In  Figure  2,  we  see  the  variation  of  {TU|  for  aodel  6  plotted  for 
six  angles  of  incidence  in  the  frequency  range  froa  0*02  to  2*0  ops* 

This  particular  aodel  was  chosen  because  it  illustrates  features  of 
the  curves  idiioh  are  found  in  all  of  the  models,  and  in  addition,  the 
frequency  variation  is  slow  enough  that  the  features  nay  be  readily 
observed. 

In  all  of  the  figures  of  this  type,  the  variation  with  frequency 
is  obvious*  With  few  exceptions,  it  has  not  been  possible  to  interpret 
the  frequency  variations  of  these  curves  as  being  due  to  any  particular 
physical  phenomenon  (for  exa^le,  multiply  reflected  waves  of  a  specif¬ 
ic  type  within  the  whole  crustal  model  or  any  seotion)*  In  this  figure, 
we  oan  observe  the  effect  of  the  angle  of  inoidenoe  on  the  frequency 
variation*  We  oan  make  the  following  observational 

1)  As  the  angle  of  incidence  increases,  the  average  amplitude 
increases*  This  might  be  expected  in  view  of  the  faot  that  the 
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dilatational  Input  will  have  a  greater  hor  is  octal  component  et  the  lar¬ 
ger  angles  of  incidence. 

2)  The  frequency  variation  is  more  rapid  at  snail  angles  of  in¬ 
cidence  than  at  larger  angles.  This  is  due  to  the  fact  that  for  the 
snail  angles  of  incidence  in  this  range,  the  P  to  S  conversion  has  a 
large  horisontal  coup  orient  which  effectively  doubles  the  niufeer  of  wave 
types  causing  interference  in  the  physical  reals. 

3)  A e  the  angle  of  incidence  increases,  the  character  of  the 
variation  changes  considerably.  Far  example,  the  peak  ihioh  occurs  at 
about  O.lU  ope  in  the  graph  for  the  29  degree  angle  of  incidence  has 
almost  disappeared  at  U9  degrees.  On  the  other  hand,  some  relatively 
large  amplitudes  at  large  angles  were  relatively  id  nor  at  s nailer 
angles.  However,  the  extremes,  that  is  the  largest  naxina  or  smallest 
mini na, tend  to  persist  and  this  suggests  their  dependence  on  the  dila- 
tational  energy. 

U)  As  the  angle  of  incidence  increases,  we  can  observe  a  system¬ 
atic  shift.  At  low  frequencies,  this  shift  is  not  so  apparent,  but  at 
the  higher  frequencies,  it  becomes  progressively  sore  evident.  The 
rate  of  shift  appears  to  increase  continuously  with  the  increase  in 
angle  of  incidence,  if  allowances  for  the  change  in  character  are  nade. 
It  will  be  seen  in  later  figures  that  a  decrease  in  the  velooity  also 
oan  cause  this  shift  to  take  place. 

In  Figure  3*  we  see  the  variation  of  |w|  for  nodal  6.  In  this 
figure,  we  oan  sake  the  following  observations  t 
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1)  As  the  angle  of  incldenoe  Increases,  the  average  aaplitude 
decreases.  This  again  would  be  expeoted  on  the  basis  of  the  vertical 
component  of  the  dilatation  being  smaller  at  large  angles  of  incidence* 

2)  The  variation  with  frequency  is  smaller  at  large  angles  of 
inoidenee  than  it  is  at  small  angles.  This  corresponds  to  the  fact 
that  for  large  anglee  of  incidence  the  vertical  component  of  P  la 
small,  and  further,  the  angle  of  inoidenee  of  the  converted  S  notion 

is  suoh  that  the  S  notion  has  only  a  snail  vertical  component  even  when 
the  angle  of  incidence  of  the  dilatational  energy  at  the  base  of  the 
crust  is  53  degrees. 

3)  As  before,  the  character  of  the  variation  changes  considerably 
when  the  angle  of  incldenoe  increases.  For  example,  the  trough  which 
appears  in  the  range  0.11*  to  0.16  ops  at  m  angle  of  29  degrees  has 
beocae  a  peak  at  k9  degrees . 

U)  Again  as  the  angle  of  incidence  increases  we  can  observe  a 
shifting  of  the  curves.  This  shift  increases  with  both  angle  and  fre¬ 
quency. 

In  comparing  Figures  2  and  3,  we  note  that  for  snail  angles  the 
average  solitudes  in  Figure  2  are  smaller  than  those  in  Figure  3.  At 
larger  angles,  the  two  are  nearly  equal.  The  frequency  variation  of 
the  |TU  |  is  more  irregular  than  that  of  |iv|.  The  shifting  effeot  is 
of  the  same  magnitude  in  both  figures. 

In  Figure  U,  we  see  the  variation  of  the  phase  of  TU  for  nodal  6, 
for  a  variety  of  anglee  in  the  same  frequency  range  aa  above.  For  this 
figure,  we  oan  make  the  following  observations! 
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1)  At  very  low  frequencies,  the  phase  shift  approaches  180  de¬ 
grees,  This  is  the  result  of  using  the  absolute  value  of  the  amplitude 
in  the  ratio  rather  than  the  horisontal  component  which  would  have  a 
sign  connected  with  it.  If  we  allow  for  this  ohange  in  phase  caused  by 
the  mathematics,  we  can  interpret  this  phase  shift  as  saying  that  at  low 
frequencies  the  incoming  wave  does  not  "see"  the  orust, 

2)  The  phase  shift  is  not  linear,  and  the  non-linearity  is  es¬ 
pecially  severe  at  small  angles  of  incidence,  where  it  often  amounts  to 
variations  of  up  to  UO  degrees.  At  larger  angles,  the  phase  shifts 
become  very  nearly  linear  in  certain  regions.  As  the  frequency  in¬ 
creases  it  appears  that  the  phase  shift  becomes  more  irregular.  How¬ 
ever,  as  we  will  see  in  later  figures,  this  is  not  progressive,  but  is 
a  periodic  phenomenon.  The  decrease  of  the  non-linearity  with  increas¬ 
ing  mgle  of  incidence  and  the  appearanoe  of  the  vertieal  phase  shifts, 
which  will  be  considered  next,  both  indicate  that  the  non-linearity  is 
due  to  the  S  conversions .  However,  we  have  been  unable  to  interpret 
this  on  the  basis  of  a  particular  ray  arrival  or  the  interference  of 
particular  ray  events, 

3)  As  the  angle  of  incidence  increases,  the  absolute  value  of 
the  slope  of  the  ourves  becomes  smaller.  This  is  related  to  the  shift 
which  appeared  in  the  amplitudes  of  TU  and  TW  .  We  can  interpret 
this  decrease  in  the  absolute  value  of  the  slope  as  the  angle  of  in¬ 
cidence  increases  in  terms  of  the  travel  tine  of  the  wave  front,  TU 
and  TW  represent  the  relations  between  two  points,  one  at  the  base  of 
the  orust,  and  the  other  at  the  free  surface  directly  above  the  first 
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point*  As  the  angle  of  incidence  increases,  the  distance  between  the 
two  points  along  the  nornal  to  the  wave  front  decreases.  If  we  de¬ 
crease  the  distance,  then  the  time  required  for  the  wave  front  to  move 
through  these  two  points  decreases,  and  the  phaae  difference  between 
two  frequencies  differing  by  an  amount  A  p  decreases  • 

Phase  difference  *  (Ap)t 
where  t  -  travel  time* 

This  explanation  was  suggested  from  a  consideration  of  plane  waves] 
and  in  the  steady-state  ease,  the  surface  cf  constant  phase  may  be  a 
curve  due  to  the  superposition  of  the  many  reflections  and  conversions 
possible*  However,  since  this  same  argument  (i.e.,  based  on  travel 
times)  explains  the  shifts  in  the  amplitude  variation  and  the  changes 
in  the  slope  of  the  phase  variation  observed  in  the  various  models,  we 
believe  that  this  explanation  is  essentially  correct* 

In  Figure  5,  we  have  the  variation  of  the  phase  of  TV  for  model  6 
for  a  series  of  angles  in  the  same  frequency  range  as  above* 

1)  At  low  frequencies,  the  phase  shift  approaches  360  degrees* 
Again  allowing  for  the  sign  convention  and  the  sign  of  the  total  ampli¬ 
tude  at  depth,  this  360  degree  phase  shift  at  low  frequencies  indicates 
that  the  low  frequencies  are  not  affected  by  the  crust* 

2)  The  phase  shift  far  TV  is  very  nearly  linear,  and  the  extent 
of  the  non-linearity  dose  not  increase  appreciably  as  the  angle  of 
inoidanoe  increases  in  this  range* 

3)  As  the  angle  of  incidence  increases,  the  absolute  value  of 
the  slope  beoomss  smaller*  This  is  due  to  the  sane  oauses  as  were 
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previously  mentioned  for  the  phase  of  TU. 

In  Figures  6,  7,  end  8,  ve  have  the  variation  in  the  frequency  de¬ 
pendence  of  |T0 1  as  a  function  of  the  model*  In  Figure  6,  the  varia¬ 
tion  is  plotted  for  an  angle  of  incidence  of  21  degrees,  Figure  7  for 
33  degrees,  and  Figure  8  for  kS  degrees*  The  observations  idiioh  will 
be  made  apply  to  all  three  figures* 

1)  The  effect  of  the  low-velocity  surfaoe  layer  in  model  1  is  to 
introduce  a  "long  period,"  large  amplitude  oscillation  in  the  frequency 
variation  of  |TU|  for  this  model. 

Further,  the  low  velooity  in  this  surface  layer  oauses  a  shift  to 
the  left  in  model  1  with  respect  to  model  2  (#xioh  is  identical  ex¬ 
cept  for  this  layer)*  This  shift  oan  be  observed  ait  frequencies  greater 
than  1*5  ope* 

2)  Vera  it  not  for  the  long  period  variation  and  the  shift  caused 
by  the  difference  in  velocities,  models  1  and  2  would  have  essentially 
the  same  frequency  variation  since  the  character  of  the  two  is  quite 
similar* 

3)  Models  2  and  3  are  similar  in  character  at  low  frequencies, 
and  although  the  character  changes  at  higher  frequencies,  it  might 
prove  to  be  practically  i^ossible  to  distinguish  between  these  two 
models  at  frequencies  less  than  0*2  ops* 

U)  Models  3  and  U  are  similar  in  amplitude,  but  the  fine  charac¬ 
ter  is  different*  This  1s  attributable  to  the  differences  in  veloci¬ 
ties*  The  shift  between  these  two  models  is  noticeable  and  is  due  to 
the  one  kilometer  difference  in  the  thickness  of  the  total  crust* 
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5)  Model  5  has  the  a am  velocities  as  nod ala  2  and  6  although 
the  ralativa  thlcknaaaaa  of  the  layers  ara  different.  The  oharaotar 
la  actually  quite  similar  in  these  three  models •  However,  the  shift 
to  tha  left  in  nodal  5  due  to  tha  relatively  large  total  cruatal 
thickness  is  so  great  as  to  oauae  tha  amplitudes  at  a  given  fraquanoy 
to  be  considerably  different.  At  higher  frequenoiee  tha  character 
also  changes  with  respeot  to  that  of  nodal  2,  but  it  still  remains 
similar. 

6)  Model  6  is  similar  to  nodal  5  in  terms  of  tha  velocities, bub 
is  much  thinner  (20  kiloneters  in  total  thickness  in  model  6  versus 
$0  kilometers  in  model  5)«  This  causes  a  considerable  shift  to  the 
right  in  the  values  of  |TU  |  in  nodel  6. 

7)  The  difference  in  the  models  at  a  given  low  frequency  is  more 
apparent  at  large  angles  of  incidence  than  at  small  angles.  This  is 
due  to  the  change  in  character,  the  change  in  amplitude,  and  the  shift. 
On  the  other  hand,  at  the  larger  angles  the  character  of  the  curves 

is  relatively  siaple.  The  character  alone  does  not  give  a  good  basis 
for  distinguishing  between  the  nodal s,  but  the  values  of  the  co¬ 
efficients  would  have  to  be  used. 

In  Figures  9,  10,  and  11,  we  have  relations  similar  to  those  shewn 

! 

in  Figures  6,  7  and  8}  however,  in  this  case  we  are  considering  | TVf |  . 
Tha  same  comments  as  were  made  about  | TU  |  are  applicable  here.  In 
ooaparing  the  two  sets  of  curves,  we  see  that  the  positions  of  the 
maxima  and  minima  of  the  two  sets  of  curves  do  not  occur  at  tha  sans 
place,  and  the  oharaoters  of  the  two  sets  ara  entirely  different* 
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In  Figures  12  through  lit,  ve  hart  the  variation  in  the  phase  of 
TU  for  the  various  models  at  three  angles  of  inoidenoe.  In  these 
figures  we  oan  observe  the  following  points t 

1)  The  phases  of  the  first  three  models  at  low  frequencies  are 
nearly  the  same. 

2)  The  relations  between  the  absolute  values  of  the  slopes  of 
the  models  are  model  6  <  model  It  <  model  3  <  model  1  <  model  2  < 
model  5* 

3)  The  difference  in  the  slopes  becomes  more  apparent  at  the  lar¬ 
ger  angles  of  inoidenoe* 

U)  There  are  definite  regions  of  non-linearity.  The  character 
of  the  non-linearity  appears  to  be  similar  in  all  of  the  models.  In 
these  regions  there  appear  to  be  dis continuities.  However,  on  the 
basis  of  model  6  we  would  expect  the  variation  to  be  continuous  al¬ 
though  very  rapid. 

In  Figures  15-17,  we  have  the  variation  in  the  phase  of  TW  for 
the  same  models  at  the  same  angles  of  inoidenoe.  Here  we  oan  make  the 
following  observations i 

1)  There  are  slight  variations  in  the  linearity  which  become 
more  apparent  at  the  larger  angles.  However,  for  most  practical  pur¬ 
poses  we  could  consider  these  phase  shifts  to  be  linear. 

2)  At  low  frequencies  the  phases  of  the  first  three  models  are 
nearly  the  same. 

3)  The  values  of  the  slopes  are  arranged  in  the  same  order  as 
were  the  phases  of  TU.  The  slopes  of  the  phase  shifts  of  TV  are  very 


figure  12.  the  effect  of  the  crustal  model  on  the  frequency  variation 
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FIGURE  16.  THE  EFFECT  OF  THE  CRUSTAL  MODEL  ON  THE  FREQUENCY  VARIATION 
OF  THE  PHASE  OF  TV  AT  AN  AN(H£  OF  INCIDENCE  OF  33  DECREES 
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FIGURE  17.  THE  EFFECT  OF  THE  CRUSTAL  MODEL  ON  THE  FREQUENCY  VARIATION 
OF  THE  PHASE  OF  TW  AT  AN  ANGLE  OF  INCIDENCE  OF  45  DEGREES 
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similar  to  the  slopes  of  the  corresponding  phase  shifts  of  TU. 

In  Figures  18-21  we  see  the  variation  of  |TU  |  ,  [TW  |  ,  the  phase 
of  TU  and  the  phase  of  TW,  respectively,  for  all  of  the  models  for  m 
angle  of  incidence  of  33  degrees  in  the  frequency  range  2.02  to  luOO 
ops.  Upon  comparing  these  figures  with  the  previous  figures,  we  see 
that  the  character  is  essentially  the  sans  as  the  corresponding  fig¬ 
ures  in  the  lower  frequency  range. 

One  noiioeable  exception  to  this  statement  1s  especially  prominent 
in  Figure  20  in  the  frequency  range  from  2.02  to  2.5  ops  for  model  1. 
Here  there  is  a  region  in  ihich  the  phase  ohanges  very  rapidly.  This 
condition  dies  out  but  does  repeat  itself  at  higher  frequencies.  We 
see  that  the  character  of  the  amplitudes  is  also  sovwhat  different 
in  this  range  (Figures  7,  10)  j  however,  this  is  not  so  noticeable.  This 
condition  is  apparently  related  to  rotational  interference  since  the 
unusual  character  appears  mainly  in  the  horisontal  component. 

Also  we  notice  that  the  shift  has  increased  to  the  point  where  it 
is  no  longer  possible  to  correlate  peaks  and  troughs.  The  values  of 
TU  and  TW  nay  be  considerably  different  at  a  given  frequency  in  this 
higher  frequency  range.  From  this  we  can  conclude  that  our  best  oppor¬ 
tunity  for  distinguishing  between  models  would  be  to  know  TU/IW  for  the 
range  of  frequency  from  0  to  3*0  ops  at  intervals  of  .015  ops.  This 
would  allow  us  to  distinguish  between  models  as  similar  as  the  first 
four  in  this  paper. 

In  order  to  test  the  validity  of  the  computations  and  to  provide 
a  comparison  with  conventional  ray  theory  in  the  real  time  domain,  the 


FIGURE  21.  THE  FREQUENCY  VARIATION  OF  THE  PHASE  OF  TV  IN  THE  FREQUENCY  RAN® 
FROM  2.02  TO  4.00  CPS  AT  AN  ANGLE  OF  INCIDENCE  OF  33  DEGREES 
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values  of  TU  and  TW  were  used  to  Fourier  synthesise  the  surfaee  wot Ion 
due  to  an  Input  pulse.  Since  the  values  of  TU  and  TV  were  known  In  the 
range  from  *02  to  10*0  ops,  an  Incident  wave  of  the  forn 

/2Q|| 

.OliTT00®  pt  dp 

was  chosen.  Upon  evaluating  the  Integral  we  have, 

INCIDENT  WAVE  -  (1/lT)  ((sin  20 rr  t)/t  -  (sin  .Ol*tT  t)/t). 

At  t-0,  this  function  has  the  vaLue  of  19*96  and  on  either  side  of  t»0 
it  appears  as  a  rapidly  decaying  sine  wave*  Because  the  first  tern  is 
Mich  larger  than  the  second  (by  a  factor  of  500),  for  small  values  of  t, 
the  predominant  features  of  this  input  wave  very  closely  resemble 
(l/TT)(8in  20 tt  t)/t.  The  value  of  this  function  does  not  ezoeed  12 
per  cent  of  the  peak  value  for  t  >  .13  seconds.  Ve  combined  this 
Input  wave  form  with  TU  and  TW  according  to  equations  (9)  in  section  2* 
These  integrals  were  evaluated  for  both  u^t)  and  wQ(t)  for  all  of  the 
models  at  all  of  the  angles.  The  time  range  which  was  chosen  for  eaoh 
angle-model  combination  was  centered  on  the  time  of  the  first  arrival 
of  the  surface  motion  as  predicted  from  independent  ray  theory  calcula¬ 
tions  which  we  carried  out.  We  choee  the  time  interval  of  the  synthesis 
at  .02  seconds.  The  synthesised  ground  motions  for  the  angles  of  25, 

33,  Ul  and  U9  degrees  are  shown  in  Figures  22-29* 

In  addition  to  confuting  the  ground  motion  from  TU  and  TV,  we 
used  the  standard  ray  theory  as  given  in  Chapters  Two  and  Three  of 
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filaetic  Wave  in  Layered  Media  (Iwing,  Jardetsky,  and  Press,  1957)  to 
oonpute  the  arrival  tins  and  amplitude  of  the  first  arrival.  In  the 
amplitude  calculations  we  took  into  aooount  the  free  surface.  In  doing 
so,  we  assumed  that  Poisson's  ratio  was  0.25  in  the  upper  layer  for  com¬ 
putational  simplicity.  In  these  models,  the  maxinun  deviation  from  this 
value  is  .01$  or  6  per  oent.  The  relative  times  and  amplitudes  of  the 
ray  theory  arrivals  are  shown  in  the  figures  by  a  solid  line.  We  see 
that  the  agreement  between  the  synthesised  motion  and  the  amplitudes  of 
the  first  arrival  based  on  ray  theory  is  very  good.  The  times  agree  to 
within  .007  seconds  and  the  amplitudes  to  within  6  per  oent.  The  author 
oonsiders  this  agreement  to  be  a  substantial  check  on  this  work.  This 
is  especially  important  since  most  of  the  calculations  cannot  be  checked 
direotly. 

In  considering  the  synthesised  motion  we  see  that  the  motion  in 
models  2-6  resembles  the  time  history  of  the  source  very  closely,  while 
the  motion  after  the  peak  arrival  in  model  1  is  quite  different.  We 
would  expect  this  because  of  the  presence  of  the  thin  surfaoe  layer  in 
model  1  whioh  is  not  present  in  the  other  models  •  The  arrival  times  of 
these  later  os o ill at ions  agree  with  those  predicted  for  P  to  S  conver¬ 
sions  at  the  uppermost  layer  and  with  multiply  reflected  P  and  S  waves 
within  this  thin  layer.  In  addition  to  the  work  shown  in  these  figures, 
we  oarried  out  more  extensive  synthesis  at  a  much  coarser  time  interval, 
and  we  found  that  the  synthesised  record  showed  energy  arriving  at  the 
arrival  times  predicted  by  ray  theory  for  multiply  reflected  and  con¬ 
verted  waves  for  all  of  the  models •  This  was  taken  as  a  further  eheok 
on  the  validity  of  these  oaloulations. 
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iu  Smmvj  and  Conclusions 

In  the  preceding  sections ,  we  here  outlined  the  use  of  the  Haskmll- 
T hoes on  matrix  method  to  determine  the  effect  of  the  crust  on  Steady- 
state  dilatational  waves.  We  showed  how  we  could  combine  a  knowledge 
of  the  transmission  effect  of  the  crust  with  the  recorded  surfaoe  motion 
to  determine  the  motion  at  the  base  of  the  orust.  Or  alternatively,  we 
showed  how  a  knowledge  of  the  surfaoe  notion  at  a  single  station  might 
be  used  to  determine  the  orustal  structure  beneath  the  station. 

Since  little  is  known  about  the  frequency  variation  of  the  trans¬ 
mission  effect  of  the  orust,  we  calculated  the  transmission  coefficients 
for  six  representative  crustal  models  for  nine  angles  of  incidence  in 
the  range  from  21  degrees  to  53  degrees.  The  frequency  variation  was 
considered  in  the  range  from  .02  to  10.0  ope.  Once  we  determined  these 
coefficients  we  used  them  to  synthesise  the  surfaoe  motion  idiich  would 
arise  from  a  pulse  of  the  form  of  (l/H )  (sin  20rrt)/t,  and  thus  obtained 
real  time  results  which  could  be  ohecked  by  o o^aris on  with  the  first 
arrivals  predicted  by  ray  theory.  From  these  computet  ions,  we  can  make 
the  following  conclusions! 

1)  Variations  in  the  total  thickness  of  the  orust  as  small  as  one 
kilometer  cause  a  noticeable  shift  in  the  positions  of  the  peala  and 
troughs  of  the  transmission  coefficients.  The  thicker  the  czust,  the 
more  rapid  is  the  frequency  variation. 

2)  For  orustal  models  having  the  same  total  thickness,  the  main 
difference  in  the  frequency  variation  of  the  transmission  coefficients 
is  in  the  character  of  the  oscillations.  At  low  frequencies  (less  than 


0*2  ops)  this  is  difficult  to  determine. 

3)  Ths  variation  in  ths  horiiontal  transmission  ooeffioimt  is 
more  irrsgular  than  that  of  ths  vertioal  coefficient  in  ths  rangs  of 
angles  of  inoidenoe  studied.  This  is  Ate  to  ths  importance  of  oocnrsrtsd 
SV  waves  in  this  range. 

U)  The  effeot  of  thin,  low  velocity  surface  layers  is  quite  pro¬ 
nounced  in  the  transnission  coefficients.  This  effeot  appears  as  a  large, 
"long  period"  oscillation  in  the  frequency  variation  of  the  traasaisalen 
coefficients. 

5)  A  law  velocity  layer  at  depth  is  very  diffioult  to  determine, 
especially  from  low  frequency  neasuresants.  However,  the  difference  la 
shift  caused  by  this  layer  is  observable,  and  if  the  transmission  co¬ 
efficients  are  known  over  a  sufficient  range  (say  .02  to  3*0  ope),  it  nay 
be  possible  to  observe  the  presence  of  the  low  velocity  layer. 

6)  In  general,  it  is  necessary  to  know  the  transmission  coeffi¬ 
cients  at  intervals  of  about  .015  ops  for  crustal  models  of  about  30 

t 

kilometers  thickness  in  order  to  observe  all  of  the  significant  features. 
For  thicker  crusts,  it  will  be  neoessary  to  know  the  values  at  even 
s seller  intervals.  Thinner  crustal  Models  could  be  investigated  at 
coarser  intervals. 

7)  The  method  of  det ernining  crustal  structure  through  the  use  of 
these  transnission  coefficients  holds  considerable  promise  when  the  mo¬ 
tion  after  the  P  arrival  shows  oosplioated  oscillations  whioh  arise  from 
the  combination  of  the  input  wavs  fora  and  the  orustal  interference. 
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8)  Whan  oosparing  the  raoorda  obtainad  from  instruments  with  vary 
different  frequency  raaponaa  ourraa  it  la  important  to  raalisa  that  tha 
aaetions  of  tha  spectrum  of  the  input  to  whioh  tha  tvo  instruments  are 
sensitive  nay  ba  very  different*  Thua  it  ia  poaaibla  in  extreme  oaaaa 
that  short  period  inatruaanta  show  substantially  larger  diaplaoaaants 
than  long  period  inatruaanta  at  tha  same  sites*  Further,  even  tha  saaa 
inatruaanta  at  diffarant  sites  at  tha  same  distance  from  tha  point  of 
generation  aay  show  different  amplitudes  due  to  the  cruatal  affect.  This 
offset  auat  ba  taken  into  account  in  such  investigations  as  source 
aeohanisa  studios  and  magnitude  determinations* 

9)  In  those  studies  in  which  frequency  dependence  is  important, 

or  in  which  the  observed  notion  is  so  aoaplicated  that  individual  events 
cannot  ba  identified,  this  asthod  is  superior  to  ray  theory  in  that  its 
application  is  more  direct.  However,  for  models  of  tha  type  studied 
here,  the  two  wthods  are  equivalent  in  that,  given  the  aouroe  fora  and 
the  crustal  nodal,  it  is  possible  to  predict  the  surface  notion  from 
either  one* 
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appendix  r 

\m  1620  COMPUTER  PROGRAM  FOR  TU  AND  TV/ 


THIS  PROGRAM  IS  DESIGNED  TO  COMPUTE  THE  TRANSFER  FUNCTIONS 
FOR  INCIDENT  L) I LATAT I ONAL  WAVES  ONLY.  . 

APPARENT  SURFACE  VELOCITY  MUST  BE  GREATER  THAN  P  VELOCITY 

in  the  top  layer  and  the  s  velocity  in  any  layer 

DIMENSION  0(60),a(60). B(60),RH0(60),DGaH(60),DGaM1(60), 
1DRA(60),uRB(60) ,0H(60) 

1030  PUNCH  1 

1  FORMAT (42HP  ORIENTED  LAYERED  SYSTEM  RESPONSE  9/25/63) 

C  READ  IN  FOUR  LETTER  IDENTIFICATION 

READ  1101.  T2 
1101  F0RMaT(a4) 

C  READ  IN  NO.  OF  LAYERS  AND  LaYER  CONSTANTS 
READ  2  #  NOL 

2  F0RMaT(I3) 

PUNCH  2, NOL 

C  READ  IN  LAYER  CONSTANTS 

DO  3  I  a  1 ,  NOL 

READ  500,  D(l)  ,A( I ) ,B( I ),RHO( I ) 

3  PUNCH  500,  o(i),a(i)»b(i ) ,RHO( I ) 

500  F0RMAT(4F 10.4) 

PUNCH  1620 

1620  FORMAT (2HWV,3X,4HANG I, 6X,4HFREQ,5X,7HAPP  vel,5x, ihu,8x,3hphu, 
18x,  1HW,8x,3HPH'7,6x,4HOMEG) 

N  *  NOL 
N1  -  NOL  -  1 

N2  -  NOL  -  2 

C  READ  IN  INITIAL  ANGLE,  INCREMENT,  FINAL  ANGLE  IN  DEGREES 
7001  READ  4  ,  aNGIP,  DaNGI.FaNGI 

4  FORMAT (3F 10.4) 

C  READ  IN  INITIAL  FREQUENCY, INCREMENT,  FINAL  FREQUENCY 
READ  4,  FREQl,OFREsi,FFREs* 

NOaNG  =  (FANG  I  -  aNGIP)/DANGI  +  1. 

NOF  =  (FFREu  -  FREq1)/DFREq  +1. 

AUG  I P  =  aNG  I  P  -  DaNG  I 
DO  310  IA  =  1, NOaNG 
FREQ  =  FREQl  -  DFRE* 

ANGIP  =  ANGIP  +  DANG  I . 

C  =  a(N0L)/SINF(AHGIP/57.295780) 

C  COMPUTE  REUSABLE  VARIABLES  FOR  THE  FIRST  LAYER 
FCOVA  =  C/a(D 
FCOVB  -  C/B(  1) 

FGAM  =  2./(FC0VB**2) 

FGAM1  -  FGaM  “1. 


FRA  =  SqRTF(aBSF(FC0Va**2  -  1.)) 

FRB  *  SQRTF(aBSF(FC0VB**2  -  1.)) 

FH  -  RHO( 1)*C*C 

C  COMPUTE  REUSABLE  VARIABLES  FOR  THE  REMAINING  LAYERS 
DO  1346  h  -  2,N0L 
COVa  =  C/a(M) 

COVB  =  C/B(M) 

OGaM(M)  -  2./ (COVB**2  ) 

DGaMI(M)  =  UGAM(M)  “1. 

DRA(M)  *  SQRTF(aBSF(COVA**2  -  1.)) 

ORB(M)  =  SQRTF(aBSF(C0VB**2  -  1.)  ) 

1346  L)H(M)  -  RHO(M)*C*C 
DO  310  IFR  »  1, NOF 
FREQ  =  FREQ  +  DFREQ 
V/VNO  -  6.283 1853*FREQ/C 
P  -  WVNO*D( 1)*FRA 

Q  »  WVNO*D( 1 )*FRB 
.  SINP  «  SINF(P) 

VI  *  SINP/FRA 
X  =  FRA*S I NP 
COSP  =  COSF(P) 

SINQ  -  SINF(Q) 

Y  -  SINQ/FRB 
Z  -  FRB*S I NQ 
COSQ  »  COSF(Q) 

C  COMPUTE  ELEMENTS  OF  A  MATRIX  FOR  FIRST  LAYER 

All  -  fgam*cosp  -  fgami*cosq 

A 12  ■  FGaM1*W  +  FGAM*Z 

A21  -  FGAM*X  +  fgami*y 

A22  »  -FGAM1*C0SP  +  FGAM*COSQ 

A31  =  FH*FGaM*FGaM1*(C0SP  -cosq) 

a32  =  FH*(FGAM1*FGAM1*W  +  fgam*fgam*z) 

a4i  =  -FH*(FGAM*FGAit*X  +  FGAM1*FGAM1*Y) 

a42  =  a31 

I F(N2)  1347,1349,1347 

c  COMPUTE  ELEMENTS  OF  A  MATRIX  FOR  REMAINING  LAYERS 

1347  UO  1345  M  =  2,N1 
GaM  =  DGaH(M) 

GaMMI  -  DGaMI(M) 
ra  =  ura(m) 

RB  -  URB(M) 

H  =  UH(M) 

P  =  wvno*d(m)*ra 

0=  WVN0*U(M)*R3 
IF(C-a(M)) 121. 122,  123 
123  SINP  =  SINF(P) 
v7  =  SINP/RA 
X  =  Ra*S I NP 
COSP  -  COSF(P) 

GO  TO  124 
122  W  -  V/VNO*u(M) 


X  -  0. 

COSP  -  1. 

GO  TO  124 
121  EXP P  -  EXPF(P) 

EXPM  -  EXPF(-P) 

V7  -  (EXPP-EXPH)/(2.*RA) 

X»  RA*(EXPM-EXPP)/2. 

COSP  »  (EXPP+EXPM)/2. 

124  S I Nq  -  SlNF(y) 

Y  -  SINq/RB 
Z  -  RB*SINy 
COSQ  =  COSF(Q) 

B 1 1  ■  GaM*COSP  -  GAMH1*C0SQ 

B 12  -  GAMM1*W  +  GAM*Z 

B 13  -  -(COSP-COS0)/H 

B14  ■  (W  +  Z)/H 

B21  -  GAH*X  +  GAMM1*Y 

B22  »  -GaHMI*COSP  +  GaM*COSQ 

B23  -  -(  X+Y)/H 

B24  *  B13 

B3i  -  h*gam*gammi*(cosp-cosq) 

B32  -  H*(GAMM1*GAMM1*W  +  GAM*GAM*Z) 

B33  =*  B22 
B34  -  B 12 

B41  -  -H*(GaM*GaM*X  +  GAHM1*GAMK1*Y) 

B42  -  B31 
843  -  B21 
B44  »  B 1 1 

multiply  matrices  ,  , 

EA 1 1  -  B 1 1*A 1 1  +  B  1 2*A  2 1  +  B13*A31  +  B14*a41 

EA 12  -  B 1 1*A 12  +  312*a22  +  B 13*A32  +  B14*a42 

EA21  =  82 1*A 1 1  +  B22*a2 1  +B23*a31  +  B24*a41 
EA22  »  B2 1*A 12  +  B22*a22  +  B23*a32  +  B24*a42 

EA31  -  B3 1*A 1 1  +  B32*a21  +  B33*a31  +  B34*a41 

EA32  *  B3 1*A 12  +  B32*A22  +B33*a32  +  B34*a42 
EA41  -  b41*a11  +  342*a21  +  B43*A31  +  B44*A4l 

Ea42  *  B41*A 12  +  942*a22  +  843*a32  +  B44*a42 

All  -  EaI 1 
Al2  =  EaI 2 
A21  -  EA21 
A22  *  Ea22 
a3  1  »  Ea3 1 
A32  =  Ea32 
a41  *  Ea41 
1345  a42  =  Ea42 
1349  a21  -  -a21 
a4 1  -  ~a41 
GaM  =  DGAM(N) 

GaMMI  ■  DGaMI(N) 

ra  -  ura(n) 
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RB  -  UR3(N) 

H  -  DH(N) 

C  COMPUTE  CLEMENTS  OF  E  INVERSE  FOR  THE  LAST  LAYER 
B 1 1  -  -  2.*(B(N)*B(N))  /(A(N)*A(N)) 

B 13  -  1./(RHO(N)*A(N)*A(N  ) 

B22  -  C*C*GAMM1/(A(N)*A(N)*RA) 

B24  -  B13/RA 

B31  =  GAMM1/(GAM*RB) 

833  -  -1./(H*GAM*RB) 

842  ■  1 . 

844  -  1./(H*GAM) 

EA 1 1  -  B11*A11  +  B 13*A3 1 
EA 12  -  B  1 1 *A 1 2  +  B13*A32 
Ea2 1  -  B22*a21  +  B24*a41 
Ea22  *  B22*a22  +  B24*a42 
Ea3 1  «  B31*Al I  +B33*A31 
Ea32  -  B31*A12  +  B33*A32 
Ea41  *  b42*a21  +  b44*a41 
Ea42  =  b42*a22  +  b44*a42 
OMEG  =  FREQ*6.2831853 

DR  -  EA21*EA32  -  Ea11*EA42  -  £Al2*EA4l  +  EA22*EA31 
Dl  -  EA1 1*Ea32+EA21*EA42-Ea12*EA31-EA22*EA41 
DENSQ  »  DR*DR  +  0 1 *D I 
UPNR  -  EA32*DI-Ea42*DR 
UPNI  *  EA32*DR  +  EA42*DI 

UPDC  -  (2./DENSQ)*SQRTF(UPNR*UPNR  +  UPNI*UPNI) 

C  UDTP  -  AMP  OF  TU, PHUPD  -  PHASE,  WDTP  »  AMP  OF  TW,  PHWPD  -  PHASE 
PHUPD  -  ATANF (UPNI /UPNR) 

PHUPD  ■  CORANG (PHUPD  , UPNI, UPNR) 

WPNI  -  EA4l*DR  +EA31*DI 
WPNR  «  -EA31*DR  +  EA4l*Dl 

WPDC  =  ( 2 .  /DENSQ) *SQRTF  (WPNR*WPNR+WPN  I  *V/PN  I ) 

PHWPD  -  ATANF (WPN I /WPNR) 

PHWPD  *  CORANG (PHWPD  ,WPNI,WPNR) 

UDTP  -  COVA*UPDC 

WDTP  -  COVA*WPOC  ,  % 

303  FORMAT (a4,f8. 3  ,F 10.5,F 10.3,F 10.4,F 10.5.F 10.4,F 10.5,f8.3) 

310  PUNCH  303, T2,aNGIP,FREQ,C, UDTP, PHUPD  ,WDTP,PHWPD  ,OMEG 
C  CONTROL  CARD  TO  CONTROL  RECYCLE 

1020  READ  7000,  CNTRL 
7000  FORMAT (F 10.2) 

I F (CNTRL )  7001, 1030,  1021 

1021  PRINT  990 

990  format ( i ihend  of  case) 

Y  =  LINK(EA31) 
end 

c  CORANG  is  a  subroutine  to  determine  the  right  quadrant 

C  FOR  THE  ANGLE  DETERMINED  BY  ATANF 

function  corang(fua,fui,fur) 

IF(FUa)  1,1,4 
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1  IF(FUI)  2,2,3 

2  ANGMUL  -  0. 

GO  TO  7 

3  ANGMUL  -  I. 

GO  TO  7 

4  IF(FUI)  5,5,6 

5  ANGMUL  —1. 

GO  TO  7 

6  ANGMUL  -  0. 

7  COUANG  -  FUA  +  ANGMUL*3. 14159265 
RETURN 

END 


} 


I 
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